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Liquid xenon is an important detection medium in direct dark matter experiments, which search 
for low-energy nuclear recoils produced by the elastic scattering of WIMPs with quarks. The two 
existing measurements of the relative scintillation efficiency of nuclear recoils below 20 keV lead 
to inconsistent extrapolations at lower energies. This results in a different energy scale and thus 
sensitivity reach of liquid xenon dark matter detectors. We report a new measurement of the 
relative scintillation efficiency below 10 keV performed with a liquid xenon scintillation detector, 
optimized for maximum light collection. Greater than 95% of the interior surface of this detector 
was instrumented with photomultiplier tubes, giving a scintillation yield of 19.6 photoelectrons/keV 
electron equivalent for 122 keV gamma rays. We find that the relative scintillation efficiency for 
nuclear recoils of 5 keV is 0.14, staying constant around this value up to 10 keV. For higher energy 
recoils we measure a value of 0.21, consistent with previously reported data. In light of this new 
measurement, the XENONIO experiment's upper limits on spin-independent WIMP-nucleon cross 
section, which were calculated assuming a constant 0.19 relative scintillation efficiency, change from 
8.8 X 10"*'' cm^ to 9.9 x lO"*"* cm^ for WIMPs of mass 100 GeV/c^ and from 4.5 x lO"** cm^ to 
5.6 X 10""" cm^ for WIMPs of mass 30 GeV/c^ 

PACS numbers; 29.40.Mc; 78.70.-g; 95.93.-|-d; 61.25.Bi 



I. INTRODUCTION 

Numerous observations point to the existence of a 
non-luminous, non-baryonic component of our universe 
known as dark matter [H-ll]- This matter could be in 
the form of a new type of particle [1, Q whose existence 
can be naturally explained as a thermal relic of the Big 
Bang. Such a particle, which would likely be massive and 
electrically neutral, is commonly referred to as a Weakly 
Interacting Massive Particle, or WIMP. There are cur- 
rently a number of efforts, worldwide, t_o directly detect 

These are 
U], fii and 



WIMPs in terrestrial particle detectors [6|-[1C 
sensitive to nuclear recoils below 100 keV 
take advantage of the fact that WIMPs, interacting pri- 
marily with atomic nuclei, would be distinguishable from 
the predominantly electromagnetic backgrounds present 
in such experiments. 

While there are varying strategies in this direct de- 
tection search, liquid xenon (LXe) has recently come to 
the forefront of the field as a powerful detection medium 
[8l-[l0|. When particles interact in LXe, the y p roduce 
prompt scintillation photons and ionization [13| which 
can be measured simultaneously to infer the energy and 
the type of interaction. The scintillation emission spec- 
trum from LXe is singly-peaked at 178 nm with a width 
of 13 nm [l^ . The scintillation yield, defined as the num- 
ber of photons produced per unit energy, depends on the 
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identity of the particle depositing the energy [T^, [l6j . 
Since WIMPs primarily interact with atomic nuclei, in 
LXe searches the quantity of interest is the scintillation 
yield of recoiling xenon nuclei. The recent XENONIO 
measurement [J [lo| has pushed the energy threshold 
of LXe detectors down to regions where the nuclear re- 
coil scintillation yield is poorly understood. This lack of 
knowledge is the source of XENONlO's largest systematic 
uncertainty. 

Measurement of the absolute scintillation yield of any 
particle species is quite difficult, so instead relative yields 
are often reported. Monoenergetic gamma rays of 122 
kcV from ^^Co are commonly used to calibrate the elec- 
tronic recoil energy scale. The scintillation yield of nu- 
clear recoils relative to that of 122 keV gamma rays is 
known as the relative scintillation efficiency, or Cef / , and 
a new measurement of this energy-dependent value is re- 
ported herein. While the absolute scintillation yield of 
electronic recoils is not linear in energy [l^l, the 122 keV 
gamma rays provide an anchor point upon which one can 
base scintillation yields of all other energies and species. 



II. EXPERIMENTAL SETUP 

A. Neutron Beam 

The determination of Ceff requires the production of 
nuclear recoils whose energies are known independent of 
their response in the LXe. Nearly monoenergetic neu- 
trons, produced in a T(p.n)"^IIe reaction, are incident 
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EJ301 Illations, and are included in the uncertainties in Table 
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Also seen in Fig. [U a 30 cni-tliick parafBn block was 
placed along the line of sight between the tritium target 
and the EJ301 scintillator, in order to block neutrons 
from directly interacting in the EJ301. In addition to 
the paraffin block, 5 cm of Pb shielded the EJ301 from 
gammas produced in the target. The data presented here 
were accumulated in several data sets with the following 
neutron scattering angles: 48°, 62°, 70.5°, and 109.5°. 



B. Detector Design 



FIG. 1: Schematic diagram of the experimental setup. In- 
coming 1 MeV neutrons scatter in the LXe detector and are 
tagged by the EJ301 neutron scintillator, for different scat- 
tering angles of 48°, 62°, 70.5°, and 109.5°. The paraffin and 
lead are used to shield the EJ301 neutron scintillator from 
direct neutrons and gamma rays. 



upon a LXe detector. Some of the neutrons which scatter 
under an angle 6 are detected by an EJ301 organic liquid 
scintillator detector [l^ (see Fig. [IJ, capable of distin- 
guishing electronic (gamma rays) from nuclear (neutron) 
recoils via Pulse Shape Discrimination (PSD) [^[lOl. In 
this way, the energy of the recoiling Xe nucleus is known 
kinematically, and is given by the relation 



(1 -cos6'). 



(1) 



where Ej. is the recoil energy. En is the energy of the 
incoming neutron, m„ and Mm are the masses of the 
neutron and nucleus, respectively, and 6 is the scattering 
angle (the approximation is valid when Afjv 3> m„ and 

En < mnC^). 

The measurements were conducted in the neutron 
beam of the Radiological Research Accelerator Facility 
at the Columbia Nevis Laboratory, also described in a 
previous study of -Ce// PH . In the present work, 1.9 
MeV protons were incident upon a tritium target, yield- 
ing 1 MeV neutrons. The T(p,n)'^He reaction produces 
neutrons over 47r sr, however, the luminosity is peaked 
in the forward direction and the energy variation due to 
the angular spread of the 1" LXe detector cell, 60 cm 
distant from the tritium target, was less than 0.09% j2^ . 
The incident proton energy, Ep, was known to within 
0.1%. These two systematic uncertainties, coming from 
the angular dependence of En and the uncertainty in Ep^ 
are considered negligible and are not included in the cal- 
culations of section IIII Bl The dominant spread in the 
incident neutron energy comes from the thickness of the 
tritium target, and was estimated to give a 1-a spread of 
±7.8% (23|. An additional component to the spread in 



comes from the finite size of both detectors, which leads 
to an uncertainty in the true scattering angle, 6. These 
uncertainties were determined by Monte Carlo (MC) sim- 




FIG. 2: Schematic diagram of the LXe detector used for the 
Ccjf measurement. Visible are four of the six PMTs used to 
view the 1 in'^ active LXe volume. 

The LXe detector was designed to allow a zero-field 
measurement of the scintillation light with very high ef- 
ficiency, covering >95% of the interior surface with pho- 
ton detectors. The design also minimizes the amount of 
LXe outside the active volume, to reduce the background 
from interactions in this passive scintillation layer. A 
schematic of the detector design is seen in Fig. [51 The 
LXe volume is viewed by six 1" -square Hamamatsu metal 
channel R8520-06-A1 photomuhiplier tubes (PMTs), sen- 
sitive to the 178 nm Xc scintillation. Four of these 
PMTs use a new bialkali photocathode, with a quan- 
tum efficiency around 40% at room temperature [2^ . 
The PMTs, held together with a polytetrafluoroethylene 
(PTFE) frame, form a cube such that each PMT win- 
dow covers a face of the cube. Both the photocathode 
and metal body of the PMTs are held at ground poten- 
tial, with positive high voltage applied to the anodes. 
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This configuration guarantees that no residual electric 
fields exist in the LXe, whose scintillation yield strongly 
depends on the applied field [1^1 (by definition, /3e// is 
the relative light yield at zero field). The PMT assem- 
bly is mounted in a stainless steel vessel, surrounded by 
fiberglass for thermal insulation. 

Xe gas, purified by the same system used in [2l|, is 
liquefied in the vessel, cooled by a copper ring cold finger 
thermally coupled to a liquid nitrogen bath. The LXe 
level is kept above the top PMT and the temperature was 
held constant at 180K (the same T and P as in XENONIO 
with fluctuations less than 0.03%. A total of 2.5 
kg of Xe was used. 

The EJ301 scintillator is contained in an aluminum 
cylinder 3" in diameter and 3" tall, viewed by a single 
Photonis XP4312B PMT and read out with the same 
electronics as the PMTs in the LXe detector. 



C. Data Acquisition and Processing 



III. DATA ANALYSIS AND RESULTS 
A. Calibrations 

The PMTs were calibrated in situ with a pulsed blue 
LED, in order to measure and monitor the gain. The light 
from the LED produces a single p.e. spectrum, whose 
mean determines the gain of the multiplier chain. With 
a complete set of such LED calibration measurements, 
the signals obtained for all acquisitions can be converted 
to a value in number of p.e. The relationship between the 
number of collected p.e. and the total number of emitted 
photons depends on the geometrical light collection effi- 
ciency, the quantum efficiency of the photocathodes, and 
the collection efficiency between the photocathode and 
the first dynodc. Although these values arc not known 
to very high precision, they represent completely linear 
processes and hence lead to a linear relationship between 
the total number of scintillation photons and the mea- 
sured number of p.e. Comparing the p.e. yields of various 
sources thus gives a measure of their relative scintillation 
yields. 



The PMT signals are fed into an amplifier, with two 
identical outputs per intput channel. One output is dig- 
itized by a 100 MHz fiash ADC, while the other output 
is fed to the triggering system. 

For the trigger, the six LXe PMT channels are com- 
bined in pairs to form three trigger channels, each dis- 
criminated at a level of 0.3 photoelectrons (p.e.). The 
logical outputs of the three discriminator channels arc 
passed to an = 3 coincidence unit. The efficiency of 
this trigger condition, determined by the Monte Carlo 
method, is -100% at 20 p.e., slowly rolHng off to ~90% 
at 10 p.e. The EJ301 trigger is taken simply as the out- 
put of the discriminated signal. 

For the measurement of the neutrons' Time of Flight 
(ToF) the LXe trigger is fed directly to the "start" in- 
put of a Time-to- Amplitude Converter (TAC), while the 
EJ301 trigger provides the "stop" after appropriate de- 
lay. The output of the TAC is digitized by the same flash 
ADC unit. Calibration of the ToF signal is discussed fur- 
ther in section UlI Al 

The shape of the signal in the EJ301 depends on the 
incoming particle species, and can be used to distinguish 
neutrons from gamma rays since the characteristic scintil- 
lation decay time is different for these particles. As a re- 
sult, the tails of pulses resulting from nuclear recoils will 
be characteristically longer than those from electronic re- 
coils. A description of the mechanisms involved in this 
process can be found in Ref. [l^. In EJ301, the "slow" 
component is two orders of magnitude longer than the 
"fast" component, reported to be 3.2 ns [l^. A PSD 
parameter is constructed by dividing the area under the 
pulse's tail by the total area of the pulse, with the tail 
defined as the part of the trace starting 30 ns after the 
peak until the trace reaches 5% of the peak value. 
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FIG. 3: (color online). The scintillation light spectrum of 122 
keV gamma rays from ^^Co, used to calibrate the electronic 
recoil energy scale. This calibration gives a scintillation yield 
of 19.64 p.e./keV. 

As Cef f is defined relative to the scintillation yield of 
122 keV gamma rays, data from a 100 fiCi ^^Co source 
were taken periodically during the experiment. Fig. [3] 
shows the spectrum from one such calibration. The ^^Co 
yield was measured to be 19.64 ± 0.07 (stat) ± 0.11 (sys) 
p.e./keVee, where the statistical uncertainty is the com- 
bination of the parameter uncertainties of the fits from 
the various calibration data, and the systematic uncer- 
tainty is taken from the variation in this yield over the 
two-day duration of the experiment. One set of PMT 
gain values was applied to all data, and thus the system- 
atic uncertainty in the ^^Co yield quoted above accounts 
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for both variations in yield and PMT gain. 

In addition to ^^Co, data were also collected from a 
^^Na source. This source emits a /?+ which promptly 
loses energy in the Na and annihilates, producing two 
511 keV gamma rays emitted simultaneously in opposite 
directions. With the source placed between the LXe de- 
tector and the EJ301 detector, the two gamma rays will 
interact at essentially the same time in the two detec- 
tors. In this way, ^^Na provides a baseline ToF=0 which, 
when used in conjunction with a variable delay generator, 
is used to calibrate the ToF measurement system. 



B. Event Selection, Backgrounds and Results 

The processing of the data acquired at each angle 
yields two parameters which can be used to select events 
of interest: the event ToF, and the PSD parameter from 
the EJ301 neutron detector. Fig. |4] shows the distribu- 
tion of events in PSD parameter and ToF. Clearly visible 
are the nuclear recoil and electronic recoil bands, in addi- 
tion to the peaks from both gamma and neutron scatters. 
The PSD cut is chosen to accept a majority of the nuclear 
recoil band while rejecting electronic recoils. The width 
of the ToF cut is 10 ns, which is the expectation based 
on the spread in £"„ and the finite size of the detectors. 
The tail of the ToF peak is due mainly to events where 
the neutron scattered in one of the detector materials 
in addition to the LXe, before interacting in the EJ301 
scintillator. Multiple scatters in the LXe also add to the 
tail, although MC simulations indicate that their overall 
contribution is less than 2%. 

Two backgrounds contribute to the LXe spectrum 
which cannot be vetoed with the cuts described above, 
and must instead be subtracted. It is clear from Fig. |4] 
that beneath the neutron peak lies a population of events 
which have a flat ToF spectrum. These are identified as 
neutrons that accidentally interacted in the EJ301 in co- 
incidence with an unrelated event in the LXe, and are 
referred to as accidentals. As these events are uniform 
in ToF space, accidentals outside of the ToF peak should 
have the same energy spectrum as those within the peak. 
The LXe spectrum of the events inside the box of Fig. 0] 
labeled "accidentals" was used as the expectation of the 
accidentals background. The region to the left of the 
peak was chosen because the peak's extended tail con- 
taminates the accidentals spectrum to the right of the 
"neutron" peak. 

The second background that cannot be vetoed comes 
from neutrons that scattered in various detector mate- 
rials in addition to the LXe, before interacting in the 
EJ301. Here referred to as materials background, MC 
simulations show that the spectrum of these events in the 
LXe follows an approximately exponential distribution in 
the region of the peak. Fig. [5ja) displays the results of 
the MC simulation of the data set at 70.5°, indicating the 
contribution from the materials background. In order to 
estimate the spectrum of these events in the real data. 
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FIG. 4: (color online). The distribution of triggered events in 
PSD vs. ToF space from the data set at 70.5°. An "upper" 
band and "lower" band are readily identifiable in the data, 
and correspond to nuclear recoils and electronic recoils, re- 
spectively. The peak at the lower left near ToF=0, due to 
gamma rays that Compton scatter in the LXe before striking 
the EJ301, is easily vetoed by the PSD cut. A population 
of accidental triggers (see text) having a flat ToF spectrum 
is visible in both bands and contributes background events 
within the neutron peak. The LXe spectra of events within 
the left box are used as the expectations of this background. 
The width of the right box — 10 ns — is chosen to accept neu- 
trons that interact in any region of the finitely-sized detectors. 



a decaying exponential was fit to the high energy por- 
tion of the distributions after subtracting the accidentals 
background. 

After applying cuts (PSD and ToF) and subtracting 
backgrounds (accidentals and materials), a spectrum re- 
sults in which the peak from single-scatter neutrons can 
be readily identified, seen as the solid circles in Fig. [51 
The horizontal scale of these spectra is given as "keVee" 
meaning "keV electron-equivalent" , indicating it is the 
energy scale derived from the ^^Co calibration. £e// is 
found from the following relation: 



-eff 



El 

En 



(2) 



where E^^ is the electron-equivalent energy (based on 
the 122 keV scintillation yield) and Enr is the true recoil 
energy. Thus, when these spectra are fit with Gaussian 
functions, the estimators of the mean, divided by the true 
recoil energy, give the Ceff values at these energies. 

The uncertainties in the recoil energies are taken di- 
rectly from the spread in the incident neutron energy 
combined with the geometrical uncertainty due to the fi- 
nite size of the detectors. These values were obtained 
from the MC simulations and are listed in the second 
column of Table HI The uncertainties in Cef / were calcu- 



5 




1 2 3 4 5 10 20 30 40 

Energy [keVee] Energy [keVr] 



FIG. 5: (color online). Selected results of the Monte Carlo simulations, which do not include the accidentals background. 
(a)-The spectrum of events tagged at 70.5°, scaled with the measured value of Cejf giving the electron-equivalent energy 
(keVee), convoluted assuming Poisson statistics for the number of p.e. and multiplied by the simulated trigger efficiency curve. 
The green histogram is the total spectrum, and the black circles indicate the true materials background. The red dashed line is 
an exponential fit to the high-energy region of the green histogram; its agreement with the true materials background confirms 
the validity of this technique's use in the real data. The shaded blue area shows the spectrum of true elastic single scattered 
neutrons. (b)-The spectrum of events tagged at 109.5°. The data are shown in the original, recoil equivalent energy scale 
(keVr) without Poisson convolution. The materials background in this region departs from the exponential behavior seen at 
lower energies, and distorts the position of the peak from true single scatters, at 20 keV. The red dashed lines are the result 
of an exponential-f Gaussian fit. The Gaussian component, centered at 22.94 ± 4.34 keV, is used as the 'true' energy of the 
Gaussian component in the real spectrum. 



lated by considering the spread in E,. mentioned above, 
statistical errors in the Gaussian fits, the variation in 
^''Co hght yield, the uncertainty in the background esti- 
mations, and the effect of the trigger threshold roll-off. 
This last uncertainty was calculated by finding the peak 
positions before and after dividing the spectra by the 
trigger efhciency discussed in section lll CI However, only 
the lowest angle (48°) was affected by this trigger roll- 
off. The asymmetric error bars of the 5 keV data point is 
due to both the trigger roll off and the actual parameter 
uncertainty in the Gaussian fit. For all angles, the dom- 
inant contribution to the uncertainty in Cef / is from the 
spread in E^. 



TABLE I: The values of £e// obtained at the four angles 
used in this study. Error bars on the recoil energies are the 
spread of En as mentioned in section [II Al combined with the 
geometrical uncertainties. The uncertainties in Ceff are the 
combination of all statistical and systematic errors mentioned 
in the text. 



e 


Er (kcV) 




48° 


5 ±0.68 


U.i'ii„o.037 


62° 


8 ±0.91 


0.137 ±0.016 


70.5° 


10 ± 1.06 


0.140 ±0.016 


109.5° 


22.94 ± 4.34 


0.205 ± 0.039 



Though the purpose of this study was to investigate 
the behavior of £e// below 10 keV, it was necessary to 
collect data from higher-energy recoils in order to estab- 
lish a connection with previous studies. For this, the 
EJ301 was placed at a scattering angle of 109.5°, corre- 
sponding to 20.0 kcV recoils. However, this angle is close 
to the minimum in the differential scattering cross sec- 
tion of 1 Me V neutrons in Xe [1^ , and so the signal from 
"true" single scatters is well below the background. Ad- 
ditionally, the materials background in this energy range 
departs from a decaying exponential. As can be seen in 
the MC data of Fig.[5l[b), the actual "bump" in the spec- 
trum, coming primarily from neutrons which have also 
scattered in the PTFE, is actually slightly higher than 
20 keV. In order to find the true energy of the peak posi- 
tion, the same procedure used in examining the real data 
was applied here to the MC data, giving a recoil energy 
of 22.94 ± 4.34 keV. The spread in Er was taken as the 
width of the Gaussian component in the MC spectrum. 



The values obtained for Cf.f f are listed in Table HI and 
additionally shown in Fig. [7] along with the results of 
previous studies [2l|, [27l - [30 | . 
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FIG. 6: Energy spectrum measured at the four angles used in this study: (a)-48°; (b)-62°; (c)-70.5°; (d)-109.5°. The data 
are shown after subtracting both accidentals and materials backgrounds. Error bars are the combined errors of the original 
spectra, accidentals, and materials background, included in the Gaussian fits indicated by the solid curves. 



IV. DISCUSSION OF RESULTS 

The data point from the measurement at 109.5° shows 
agreement with other measurements whose high-energy 
behavior averages out to Ceff ~ 0.19. Below 10 keV, our 
values are substantially lower than the central values of 
Chepel et al. [l^j, with a considerable improvement in 
precision. The central value at 10 keV is consistent with 
the lowest-energy data point of Aprile et al. [2l| , enforc- 
ing the accuracy of this measurement. Unfortunately, the 
theoretical models of neither Lindhard [HI nor Hitachi 
[l6| can shed any light on the behavior of Cef / in this en- 
ergy range. Hitachi's model, which attempts to take into 
account incomplete charge recombination and additional 
electronic quenching, is based on Lindhard quenching as 
well as the Thomas-Fermi approximation; for Xe nuclear 
recoils, both break down below 10 keV [32l. Issj. 



As mentioned in the introduction, the uncertainty in 
Ceff at low recoil energies presents the largest system- 
atic uncertainty in the results of the XENONIO dark 
matter experiment, where it was chosen to use a flat 
Ceff = 0.19 as a compromise between the seemingly op- 
posing trends observed by Chepel and Aprile. Under this 
assumption, the WIMP-nucleon spin-independent cross 
section for WIMPs of mass 100 GeV/c^ was constrained 
to be less than 8.8 x lO"'''* cm^ and 4.5 x lO""*"^ cm^ at 
30 GeV/c^, indicated by the solid blue curve in Fig |8l 
Allowing for Cef / scenarios below 20 keV that cover the 
values allowed by both Chepel and Aprile gives upper 
limits that vary by ~40% at 30 GeV/c^ and ~18% at 
100 GeV/c^, with variations becoming less severe with 
increasing WIMP mass. With an Ceff model that fol- 
lows the new data points of this study, the resulting up- 
per limit is shown in Fig [5] as the blue dashed curve. The 
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FIG. 7; (color online). Measured £e// values as a function 
of Xe nuclear recoil energy. Symbols correspond to (o)-tliis 
work; (□)-Che2el et al. [l^l; (A)-Aprile et al. [IJ; (0)- 
Akimov et al. [2^; (x)-Bernabei et al. [29| : (v)-Arneodo et 
al. [s^. The solid gray curve is the result from a recent best- 
fit analysis of XENON 10 AmBe source data and MC [3l|. 
Also shown is the theoretical prediction of Hitachi (dashed 
line) [ill. 



limit is shifted up to 9.9 x 10 cm^ and 5.6 x 10 '^'^ 
cm^ for WIMPs of mass 100 GcV/c^ and 30 GcV/c^, 
respectively. 

It has become clear from XENONIO that future dark 
matter searches using LXe must have sensitivity to nu- 
clear recoils below 10 keV in order to be competitive. The 
improved understanding of Ceff's behavior presented in 
this study not only permits a more precise interpretation 
of XENONlO's results, but benefits future dark matter 
searches also using LXe. Several next generation LXe 
dark matter searches are currently in operation or un- 
der construction, such as XENONIOO HjI, LUX [H and 
XMASS [1^. These experiments will begin to probe for 
the first time those regions of parameter space most fa- 
vored by many theoretical models, and will consequently 
rely quite heavily on a precise understanding of LXe's 
scintillation efficiency for low energy nuclear recoils when 
interpreting their results. This is true in the case of a null 
result and especially in the case of a positive signal. If 
and when such a signal is detected, a measurement of 
the WIMP mass, for example, which relies on analyzing 
the energy spectrum of recoils, will be affected by the 
precision to which Cef / is known. 



SUMMARY 




WIMP Mass [GeV/c"] 



FIG. 8: The upper limit on the WIMP-nucleon spin- 
independent cross section based on the 58.6 live days of 
XENONlO's WIMP search, shown with a flat Ceff= 0.19 
(solid). An £e// function consistent with the results of this 
study, applied to the same XENONIO data is shown as well 
(dashed). 



direct detection with liquid xenon detectors. We show 
that at 10 keV and below, this efficiency is lower than 
the average value of Cef f =0-19 while the measurement 
in the literature |27| suggests a rise in Cef / at these en- 
ergies, albeit with large errors. In light of the results of 
our measurement, the XENONIO spin-independent limit 
is shifted up for WIMPs of mass 100 GeV/c^ by 12.5%, 
while the high-mass regime is relatively unchanged. 



The work presented here represents a new measure- 
ment of Cef f in an energy range where it is poorly under- 
stood, but highly important to the field of dark matter 
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